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In  this  report,  we  present  the  basic  analytical  MOSFET  model  which  describes  both 
the  below  and  above-threshold  regimes  of  device  operation.  The  description  is 
based  on  a  change  control  model  which  uses  one  unified  expression  for  the 
effective  differential  channel  capacitance.  The  model  also  accounts  for 
series  drain  and  source  resistances,  velocity  saturation  in  the  channel, 
finite  output  conductance  in  the  saturation  regime,  and  for  the  threshold 
voltage  shift  due  to  drain  bias  Induced  lowering  of  the  injection  barrier 
between  the  source  and  the  channel  (DIBL) .  The  model  parameters,  such  as 
the  effective  channel  mobility,  the  saturation  velocity,  the  source  and 
drain  resistances, etc.  are  extractable  from  experimental  data.  This  makes 
the  model  very  suitable  for  Inclrporation  into  our  circuit  simulator,  AIM-Spice. 

We  apply  the  characterization  procedure  based  on  this  model  to  a  MOSFET  with 
a  quarter  micron  gate  length  and  obtain  excellent  agreement  with  experimental 
data.  The  approach  developed  in  this  report  was  also  applied  for  modeling  other 
field  effect  transistors,  such  as  GaAs  Metal  Semiconductor  Field  Effect  Transistors 
and  Heterostructure -Field  Effect  Tr«n»iet€Frs. 
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I.  Abstract. 

The  objective  of  this  program  was  to  develop  a  model  fo  modeling 
radiation  effects  in  CMOS.  In  this  report,  we  describe  a  new  MOSFET  model 
suitable  for  modeling  radiation  effects  and  demonstrate  the  application  of  this 
model  to  commercial  CMOS  before  and  after  radiation  and  after  room 
temperature  annealing. 

II.  Unified  MOSFET  model. 

In  this  Report,  we  present  the  basic  analytical  MOSFET  model  which 
describes  both  the  below  and  above-threshold  regimes  of  device  operation.  The 
description  is  based  on  a  charge  control  model  which  uses  one  unified  expression 
for  the  effective  differential  channel  capacitance.  The  model  also  accounts  for 
series  drain  and  source  resistances,  velocity  saturation  in  the  channel,  finite 
output  conductance  in  the  saturation  regime,  and  for  the  thresnold  voltage  shift 
due  to  drain  bias  induced  lowering  of  the  injection  barrier  between  the  source 
and  the  channel  (DIBL).  The  model  parameters,  such  as  the  effective  channel 
mobility,  the  saturation  velocity,  the  source  and  drain  resistances,  etc.  are 
extractable  from  experimental  data.  This  makes  the  model  very  suitable  for 
incorporation  into  our  circuit  simulator,  AIM-Spice  (see  Ytterdal  et  al.  (1991)). 
It  is  also  quite  useful  for  MOSFET  design  and  characterization.  We  apply  the 
characterization  procedure  based  on  this  model  to  a  MOSFET  with  a  quarter 
micron  gate  length  and  obtain  excellent  agreement  with  experimental  data. 

The  approach  developed  in  this  Report  was  also  applied  for  modeling  other 
field  effect  transistors,  such  as  GaAs  Metal  Semiconductor  Field  Effect 
Transistors  and  Heterostructure  Field  Effect  Transistors  (see  Fjeldly  and  Shur 
(1991)). 

1.  Unified  Channel  Capacitance 

The  charge  control  in  the  MOSFETs  can  be  described  in  terms  of  a  unified 
expression  for  the  channel  capacitance,  valid  for  all  bias  voltages.  As  was 
discussed  Jeng  (1988),  the  differential  gate-channel  capacitance  per  unit  area  can 
be  written  as 
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is  the  sheet  density  of  inversion  charge  at  threshold. 

A  unified  carrier  sheet  charge  density  can  now  be  expressed  in  terms  of  the 
effective,  unified  differential  channel  capacitance.  For  an  n-channel  MOSFET  at 
very  low  drain-source  bias  (well  below  saturation),  we  find 


=  2qn^ 

Since  the  drain  current  is  assumed  to  be  small,  the  intrinsic  gate  voltage  swing, 
Vgt,  in  this  expression  can  be  replaced  by  its  extrinsic  counterpart,  Vgt  =  Vgt  + 
IdRs^  where  Id  is  the  drain  current  and  Rs  is  the  source  series  resistance. 

The  capacitance  model  described  above  can  be  used  for  calculating  the 
channel  contribution  to  the  gate-source  capacitance,  Cgs,  and  the  gate-drain 
capacitance,  Cgd,  by  using  an  approximation  similar  to  that  used  in  the  model  by 
Meyer  (1971): 
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Here,  Cgc  =  CgcLW,  Vsaie  =  hatigeh  is  the  effective  extrinsic  saturation  voltage, 
I  sat  is  the  saturation  drain  current,  gch  is  the  extrinsic  channel  conductance  at  low 
drain-source  bias,  and  Vdse  is  an  effective  extrinsic  drain-source  voltage.  V^se  is 
equal  to  Vds  for  Vds  <  V sat  and  is  equal  to  V sate  lor  V ds  >  ^ sate-  In  order  to 
obtain  a  smooth  transition  between  the  two  regimes,  we  interpolate  Vdse  by  the 
following  equation; 
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where  me  is  a  constant  determining  the  width  of  the  transition  region  between  the 
linear  and  the  saturation  regime.  The  capacitance  C/in  eqs.  (6)  to  (7)  is  the  side 
wall  and  fringing  capacitance  which  can  be  estimated,  in  terms  of  a  metal  line  of 
the  length  W,  as 


9  =  Pc  (9) 

where  Be  is  on  the  order  of  0.5  (see  Gelmont  et  al.  (1991)). 

The  source-substrate,  drain-substrate,  and  gate-substrate  capacitances  also 
affect  MOSFET  performance  in  a  circuit.  Our  model  (implemented  in  AIM- 
Spice)  uses  standard  expressions  for  these  capacitances  which  can  be  found,  for 
example,  in  Shur  (1990),  p.  393. 

2.  Current-Volfage  Characteristics 

At  drain-source  voltages  well  below  the  saturation  voltage,  the  drain 
current  can  always  (i.e.,  at  both  above  and  below  threshold)  be  expressed  as 
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where  the  intrinsic  channel  conduction  in  the  linear  region,  gchi^  is  related  to  its 
extrinsic  counterpart,  gch^  as  follows; 


(11) 


Here,  Ri  =  +  Rj  is  the  sum  of  the  source  and  the  drain  series  resistances.  (In 
deriving  eq.  (11),  we  assumed  that  the  drain  current  is  so  small  that  the 
transconductance  is  much  less  less  than  channel  conductance).  The  intrinsic 
channel  conductance  can  be  written  as 
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where  %  is  the  unified  carrier  sheet  density  at  low  drain-source  voltage  given  by 
eq.  (5)  and  [in  ic  the  low-field  carrier  mobility.  In  eq.  (10),  Vps  and  ^ds  are  the 
intrinsic  and  extrinsic  drain-source  voltages,  respectively,  which  are  related  by 
^ds  =  ^DS  +  We  note  that,  in  the  subthreshold  regime,  the  channel 

resistance  will  normally  be  much  larger  than  the  parasitic  source  and  drain 
resistances  such  that  gch-gchi>  and  Vgt-VcT- 

The  MOSFET  drain  saturation  current  for  the  above-threshold  regime  was 
derived  for  the  two-piece  velocity  saturation  model  (see,  for  example,  Shur 
(1990)).  This  current  can  be  expressed  as 
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Here,  Vl  =  FsL  where  Fs  is  the  saturation  field  and  L  is  the  gate  length. 

In  order  to  link  the  linear  and  the  saturation  regions  by  one  single 
expression,  we  propose  to  use  the  following  basic  interpolation  formula  for  our 
model: 
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where  m  is  a  parameter  that  determines  the  shape  of  the  characteristics  in  the 
knee  region. 

We  should  notice  that  this  model  can  utilize  different  approximations  for 
hat-  For  examnle,  we  can  use  the  model  proposed  by  Sodini  et  al.  (1984)  and 
utilized  in  the  popular  MOS  model  BSIM  (see  Jeng  ei  al.  (1987)).  Still  another 
alternative  is  to  use  the  empirical  model  proposed  by  Sakurai  et  al.  (1991)  {hat  = 
BW{Vg^  -  Vt)'^JL  where  B  and  n  are  empirical  fitting  parameters).  In  this  case, 
one  can  use  the  analytical  equations  for  the  propagation  delays  of  logic  gates  by 
Sakurai  et  al. 

We  now  want  to  improve  this  basic  model  to  account  for  various  non-ideal 
physical  mechanisms  such  as  subthreshold  current,  drain  bias  induced  barrier 
lowering  (DIBL),  and  other  factors  that  contribute  to  a  finite  output  conductance 
in  saturation. 

The  subthreshold  l-V  characteristics  of  the  MOSFET  in  the  long  channel 
limit,  these  characteristics  can  be  expressed  as 
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where  On  =  {J-n^th  is  the  carrier  diffusion  coefficient.  From  eq.  (15),  we  find 
the  following  subthreshold  saturation  current: 
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(The  corresponding  saturation  voltage  is  of  the  order  of  the  thermal  voltage, 

Vth-) 

A  unified  expression  for  the  saturation  current,  valid  both  above  and  below 
threshold,  can  be  obtained  by  substituting  the  gate  voltage  swing  Vgt  in  eq.  (15) 
by  an  equivalent  gate  voltage  swing  which  coincides  with  Vgt  well  above 
threshold  and  is  equal  to  2Vf/i  below  threshold: 
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Here,  5  determines  the  width  of  the  transition  region.  Typically,  5  =  3  is  a  good 
choice.  The  unified  saturation  current 

V 

j _ \hi  gtfi 

reverts  to  eq.  (13)  above  threshold  and  reduces  to  eq.  (16)  below  threshold. 
Hence,  using  eq.  (18)  for  the  saturation  current,  eq.  (14)  becomes  a  unified 
expression  for  the  l-V  characteristics  with  the  correct  limiting  behavior  in  both 
the  above  and  below  threshold  regimes. 

An  important  effect  which  must  be  taken  into  account  for  an  adequate 
description  of  the  MOSFET  behavior  is  the  dependence  of  the  threshold  voltage 
on  the  drain-source  voltage.  This  dependence  can  be  fairly  accurately  described 
by  the  equation 
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where  Vtq  is  the  threshold  voltage  at  zero  drain-source  voltage  and  o  is  a 
coefficient  which  may  depend  on  the  gate  voltage  swing.  This  translates  into  the 
following  shift  in  Vgti 
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ds 
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where  Vgto  is  the  gate  voltage  referred  to  threshold  at  zero  drain-source  bias. 
The  effects  of  the  drain  bias  induced  threshold  voltage  shift  can  be  accounted  for 
in  the  /-V  characteristics  of  eq.  (14))  by  using  eq.  (20)  in  the  expressions  for  the 
the  saturation  current  and  the  linear  channel  conductance.  This  will  result  in  a 
finite  output  conductance  in  the  saturation  region.  However,  as  discussed  in  by 
Klassen  et  al.  (1988),  a  will  depend  on  the  gate  voltage  swing  since  the  effect  of 
drain  bias  induced  barrier  lowering  (DIBL)  will  vanish  well  above  threshold. 


Hence,  for  modeling  purposes,  we  adopted  the  tollowing  empirical  expression  tor 


1  +  exp 


which  gives  a  — »  Oo  for  V gt  <  V^i  and  o  — >  0  for  V g[  >  Vat-  The  voltage 
Va  determines  the  width  of  the  transition  between  the  two  regimes. 

For  a  more  accurate  description  of  the  saturation  region,  it  is  possible  to 
include  an  additional  empirical  factor.  (1  +  in  the  expression  for  the  I-V 

characteristics  (as  usually  done  in  SPICE  FET  models).  The  parameter  X  is 
related  to  physical  effects  that  are  not  explicitly  included  in  our  model,  such  as, 
for  example,  gate  length  modulation.  It  is  important  to  emphasize  that  X  is  a 
strong  function  of  the  gate  length,  the  aspect  ratio,  the  contact  depth,  and  other 
factors.  As  a  rule  of  thumb,  we  can  assume  X  ~  kolL  where  ko  is  a  constant. 
Hence,  our  model  l-V  characteristics  can  finally  be  expressed  as 
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where 


sate  ^sat^^ch 

is  the  effective,  extrinsic  saturation  voltage  (the  same  as  used  in  eqs  (6)  -  (8)). 
When  the  experimental  l-V  characteristics  have  a  finite  output  conductance,  gc/is. 
in  saturation,  the  saturation  current  is  conveniently  defined,  for  the  purpose  of 
comparison  with  the  present  model,  as  the  the  interReport  of  the  linear 
asymptotic  behavior  at  small  drain-source  voltages  {Id  -  gchVds)  and  in  deep 
saturation  {Id  =  gchVsate  +  gchs^ds)-  Hence,  we  can  express  V^ate  in  terms  of  the 
saturation  current  as:  Vsaie  =  (1  -  gchsIgchVsatIgch-  Normally,  for  well-behaved 
transistors,  gchs  «  gch  and  we  can  use  eq.  (23).  However,  at  sufficiently  large 
gate  voltage  swing,  i.e.,  for  Vgr>  Vot  +  3^0,  the  effects  of  the  drain  voltage 
induced  shift  in  the  threshold  voltage  vanishes  and  we  obtain  gchs  *  ^gchYsate  in 
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deep  saturation.  This  can  be  used  to  deterrnine  A,  as  part  of  the  device 
characterization  and  to  obtain  the  corrected  relationship  Vsate  =  igchihat  + 
to  be  used  in  eq.  (22)  well  into  strong  inversion. 

The  quality  of  the  present  unified  MOSFET  model  is  illustrated  in  Figs.  1 
and  2  for  a  deep  submicrometer  NMOS.  Fig.  1  shows  the  above-threshold  I-V 
characteristics  and  Fig.  2  shows  the  subthreshold  characteristics  in  a  semilog 
plot.  As  can  be  seen  from  the  figures,,  our  model  quite  accurately  reproduces  the 
experimental  data  in  the  entire  range  of  bias  voltages,  over  several  decades  of  the 
current  variation.  A  residual  drain  current  at  large  negative  gate-source  voltages 
is  not  accounted  for  in  our  model.  This  current  is  caused  by  a  small  leakage  of 
the  reverse  biased  p-n  drain-substrate  junction.  The  parameters  used  in  this 
calculation  were  obtained  using  the  parameter  extraction  procedure  described 
below  (see  Section  3). 


drain-source  voltage  (V) 


Fig.  /.  Above-threshold  experimental  (clots)  and  calculated  (solid  lines)  I-V  characteristics  for  a 
deep  submicrometer  NMOS  with  effective  gate  length  L  =  0.25  pm.  Device  parameters  are: 
effective  gate  width  =  10  pm,  oxide  thickness  dj  =  5.6  nm,  mobility  pn  =  0.028  m^/Vs, 
saturation  velocity  Vj  =  4xl(f  mis,  saturation  knee  parameter  m  =  4.0,  source  and  drain  series 
resistances  Rg  =  Rd  =  75  ohm,  output  conductance  parameter  X  =  0.048  V^,  subthreshold 
ideality  factor  t\=1.32,  threshold  voltage  at  zero  drain-source  bias  0.44  V,  threshold 

voltage  parameters  Gq  =  0.048,  \(j=  0.2  V,  W(ji  =  1 .7  V.  Substrate-source  bias  Vbs  =  0  K 
Top  curve:  Vgs  =  3.0  V,  step:  ■0.5  V.  (Datafrom  Chung  et.al.  (1991).) 


{•ate-source  voltage  (V) 


Fig.  2.  Subthreshold  experimental  (  s  mhol.s)  and  calculated  (solid  lines)  I-V  characteristics  for  a 
deep  submicrometer  NMOS  with  effective  gate  length  L  =  0.25  pm  (the  same  as  shown  in  Fig. 
I).  It  denotes  the  threshold  current. 

The  advantage  of  the  model  described  in  this  Report  is  its  simplicity  and 
applicability  to  short  channel  devices.  However,  it  is  also  important  to  point  out 
the  limitations  of  this  simple  model.  For  example,  the  output  conductance  of  long 
channel  devices  in  the  saturation  regime  is  equal  to  zero  (ideal  saturation).  The 
only  way  to  approximate  such  an  ideal  saturation  using  eq.  (14)  (or  eq.  (22)  with 
A,  =  0)  is  to  choose  . a  large  value  of  the  shape  parameter  m.  However,  this  leads 
to  a  very  abrupt  transition  from  the  linear  to  the  saturation  regime  in  the 
computed  current-voltage  characteristics.  In  order  to  obtain  a  model  which 
describes  accurately  devices  with  all  channel  lengths,  we  may  choose  to  employ 
different  equations  for  the  linear  and  the  saturation  regimes  while  preserving  the 
continuity  of  the  drain  current  and  its  derivatives  with  respect  to  the  terminal 
voltages.  Such  an  approach  was  considered  by  Park  et  al.  (1991)  and  Moon  et  al. 
(1991). 

The  dependence  of  the  field  effect  mobility  on  the  gate-source  voltage  can 
be  accounted  for  within  the  framework  of  the  model  considered  in  this  Report. 
However,  the  effective  low  field  mobility  may  also  depend  on  the  drain-source 
voltage  because  the  channel  potential  affects  the  surface  electric  field  and,  hence, 
the  surface  scattering.  To  some  extent,  these  effects  can  be  indirectly  accounted 
for  in  the  present  model  by  adjusting  the  values  of  the  parameters  X  and  m. 
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However,  the  ability  of  the  model  to  reproduce  these  effects  is  limited.  Here 
again,  the  more  sophisticated  models,  such  as  the  models  developed  by  Park  et  al. 
(1991)  and  Moon  et  al.  (1991).  (These  models  have  already  been  implemented  in 
AIM-Spice.  However,  it  is  hard  to  compete  with  simplicity  and  with  a  set  of 
equations  that  is  easy  to  understand.  Besides,  for  most  practical  modem  day 
CMOS  technologies,  the  model  described  in  this  Report  is  quite  adequate. 

The  usefulness  of  the  present  unified  MOSFET  model  is  illustrated  in  the 
next  Section  where  we  describe  the  parameter  extraction  technique  using  this 
model. 

3.  Parameter  Extraction 

In  order  to  apply  the  pre.sent  model  in  circuit  simulation,  it  is  necessary  to 
find  a  practical  way  of  extracting  the  model  parameters  from  experimental  data. 
The  accuracy  of  the  parameter  extraction  depends  on  the  accuracy  of  the  model 
and  on  the  accuracy  and  completeness  of  the  experimental  data. 

In  real  life,  there  are  usually  problems  with  both  the  model  and  the 
experimental  data.  Even  though  the  analytical  model  developed  in  this 
Reportreflect  many  important  effects  occurring  in  MOSFETs,  such  as  velocity 
saturation,  parasitic  series  resistances,  the  dependence  of  the  threshold  voltage  on 
the  drain  voltag;e,  and  finite  output  conductance,  many  other  important 
phenomena  are  included  indirectly,  at  best.  For  example,  our  model  treats  the 
dependence  of  the  mobility  on  the  bias  conditions  only  through  effective  values 
and,  indirectly,  through  the  fitting  parameter  X.  Furthermore,  many  complicated 
phenomena  associated  with  the  threshold  voltage  dependence  on  the  drain-source 
voltage  are  accounted  for  in  a  set  of  simple,  empirical  formulae  (see  eqs.  (19)  to 
(21)). 

Hence,  our  approach  to  the  device  characterization  of  MOSFETs  is  to 
extract  effective  device  parameters,  such  as  the  threshold  voltage,  Vto^  the 
mobility,  \in,  the  saturation  velocity,  vy,  the  subthreshold  ideality  factor,  q,  the 
source  and  drain  series  resistances,  Rs  and  Ruy  the  output  conductance  parameter, 
the  knee  shape  parameter  m,  the  threshold  voltage  coefficient.  Go,  and  the 
voltages  Vat  and  Va  characterizing  the  dependence  of  parameter  G  on  the  gate 
voltage  swing.  In  other  words,  we  describe  a  parameter  extraction  procedure 
which  gives  fairly  reasonable  results  without  parameter  adjustment  and  which  can 
produce  a  nearly  perfect  fit  with  some  parameter  adjustment.  We  should  notice 
that,  even  without  a  precise  knowledge  of  the  device  geometry  and  process 


parameters,  we  can  obtain  a  good  fit  to  the  device  characteristics  and.  hence, 
predict  the  circuit  performance  fairly  accurately  using  the  AIM-Spice  simulation 
program.  Moreover,  we  still  retain  a  limited  ability  to  predict  how  device 
characteristics  will  scale  with  changes  in  the  doping  level  and  profile  and  device 
geometry.  In  other  words,  our  model  and  parameter  extraction  procedure 
attempt  to  give  the  best  possible  answer  based  on  the  available  information. 

As  an  example,  we  perform  the  parameter  extraction  for  the  device  with 
the  characteristics  shown  in  Fig.  1  and  2.  We  start  by  extracting  the  drain  bias 
dependence  of  the  threshold  voltage  from  the  measurements  of  the  drain  current 
versus  gate-source  voltage  for  different  drain-source  voltages.  Plotting  the 
results  of  these  measurements  in  a  semilog  scale,  as  illustrated  in  Fig.  2,  we  can 
determine  riV'^;,  frorn  the  slope  of  the  characteristics  in  the  subthreshold  regime. 
We  then  estimate  the  threshold  current  in  saturation,  /f,  from  eq.  (16)  using  eq. 
(4)  for  electron  sheet  density  at  threshold,  i.e.. 


t  (^d.+^d)L 


(24) 


At  this  point,  we  do  not  need  to  know  accurate  values  for  the  electron  mobility, 
jin,  the  effective  channel  length,  L,  and  the  oxide  thickness,  di,  because  the 
threshold  voltage  is  fairly  insensitive  to  the  exact  value  of  I [.  (For  our  device.  It 
*  3x10'^  A.)  The  intercept  of  the  measured  curves  with  the  line  Id  =  It  yields 
the  values  of  the  threshold  voltage  versus  drain-source  voltage  (see  Fig.  2).  A 
linear  interpolation  of  this  dependence  yields  the  parameter  a<7.  For  the  present 
device,  we  find  Vto  =  0.44  V  and  Oo  ~  0.048. 

The  next  step  is  to  find  the  shape  parameter,  m,  for  the  knee  region  of  the 
characteristics.  For  near  ideal  characteristics,  with  low  output  conductance  in 
saturation,  the  extraction  method  is  as  indicated  in  Fig.  3.  Each  characteristic  is 
approximated  by  two  linear  pieces,  one  describing  correctly  the  linear  region  at 
small  drain-source  voltages  (slope  gch)  and  one  approximating  the  characteristic 
in  the  saturated  region  (slope  gchs)-  The  linear  pieces  intersect  at  the  saturation 
current,  Isat,  at  the  drain-source  voltage  voltage  Vds  =  gchUsat  *  Vsate  (see  the 
discussion  in  Section  2).  According  to  eq.  (22),  the  real  current  at  this  drain- 
source  voltage  is  Id  =  hat  -  A/  =  Isatl'2-^^’^-  Hence,  by  locating  the  interReport  in 
the  two-piece  linear  model  applied  to  the  experimental  characteristic,  as  shown  in 


Fig.  3,  and  measuring  the  currents  I  sat  and  l^at  -  we  can  determine  m  as 
follows: 


ln(  2) 

(25) 

For  the  curve  shown  in  Fig.  3,  we  obtain  m  =4. 


Fig.  3.  Method  for  extracting  the  shape  parameter  m.,  the  saturation,  current,  Isat-  the 
parameter  X.  The  I-V  characteristic  is  taken  from  Fig.  I  ('Vgs  =  3.0  V). 

From  the  output  conductance  in  deep  saturation,  gchs>  for  the  largest  value  of 
Vgs,  we  also  determine  the  value  of  parameter  X  from  the  expression 

X  =  _  ^(-'hs 

^ch^ sate  ^  sat 

(see  the  discussion  in  Section  2).  From  Fig.  3,  we  find  A.  *  0.048  V'^ 

The  next  step  involves  the  determination  of  the  total  series  resistance,  Rt  = 
Rs  +  Rd*  and  the  transconductance  parameter 

e. 

£,(</.  +A</)  (27) 

In  the  above-threshold  regime,  the  intrinsic,  linear  channel  conductance  can  be 
expressed  as  follows  in  terms  of  P:  gchi  -  ^Vgt  *  ^Vgf.  Hence,  we  can  write 
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ch 


'  chi 


+  R. 


at 


(28) 


This  allows  (3  and  to  be  determined  from  the  slope  and  the  intercept  with 
^/8ch  =  0  by  plotting  l/gch  versus  Vi>[.  Such  a  plot  is  shown  in  Fig.  4  for  the 
NMOS  device  under  investigation,  resulting  in  the  following  parameter  values: 
Rf  ~  150  ohm  and  (3  =  O.Olt  A/V-.  However,  we  note  that  since  the  series 
resistance  in  this  short  channel  MOSFET  dominates  the  overall  device  resistivity 
in  the  linear  regime,  the  determination  of  (3  by  this  method  is  somewhat 
uncertain.  Besides,  this  technique  of  extracting  (3  emphasizes  the  properties  of 
the  I-V  characteristics  in  the  linear  regime  while,  in  fact,  the  properties  of 
saturation  region  are  of  greater  interest  in  digital  applications.  Therefore,  we 
return  below  with  alternative  ways  of  estimating  (3. 


Fig.  4.  Plot  of  inverse  channel  conductance  versus  inverse  gate  voltage  swing,  used  for  the 
determination  of  the  parameters  Rt  and  fi. 


For  MOSFETs,  we  normally  have  Rs  Rd  ~  Rt/'^..  If  there  is  an 
asymmetry  between  the  source  and  drain,  an  accurate  estimate  of  the  difference 
between  Rs  and  Rd  can  be  made  by  measuring  the  normal,  extrinsic  device 
transconductance  in  the  saturation  regime,  gm,  and  its  counterpart,  gmrt  with  the 
source  and  drain  interchanged  (at  constant  drain  current).  Since  the  intrinsic 
transconductance  is  the  same  in  both  cases,  we  find 

-  V" 

a  dffi  6mr 


(29) 


In  our  our  case,  however,  we  used  =  Rd  =  Rt^2  =  75  ohm. 
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Fig.  5.  Measured  saturation  current  (defined  in  Fii^  J)  versus  gate-source  voltage  for  the  same 
device  as  in  Fig.  I.  This  plot  is  used  to  determine  the  e.arinsic  transconductance  in  saturation 


Fig.  5  shows  the  dependence  of  the  above-threshold  saturation  current,  Isah  on 
the  gate-source  voltage.  From  a  suitable  interpolation  formula  representing  the 
experimental  data,  we  can  determine  the  extrinsic  transconductance  in  saturation, 
gm,  as  a  function  of  V^s  by  differentiation.  The  intrinsic  transconductance  is 
subsequently  obtained  as  gmi  =  Hm/i  I  *  }^mR.s)-  According  to  our  model,  we  can 
write 


2  ~  2  ^  2  2 

P  p  (30) 

where  Vqt -  f^shat-  Hence,  by  plotting  l/gmi^  versus  we 

obtain  P  from  the  slope  and  from  the  intercept  with  the  l/gmi^  axis  (at 

1/Vg7^  0).  A  plot  based  on  eq.  (30)  and  the  data  in  Fig.  5  is  shown  in  Fig.  6. 

As  can  be  seen  from  this  figure,  the  slope  in  this  particular  case  is  too  small  to 
serve  as  a  reliable  estimate  for  p.  However,  from  the  intercept,  we  find  that 
PV/^  »  0.026  AfV  which  corresponds  to  a  saturation  velocity  *  4x1  O'*  m/s. 
This  value  of  seems  to  be  too  low  compared  to  the  well  established  measured 
value  of  Vs  »  8x10"*  m/s  (see,  for  exaple,  Chan  et  al.  (1990)).  However,  we 
checked  that  the  values  of  extracted  from  these  experimental  data  by  other 
models,  such  as  the  model  presented  by  Park  et  al.  (1991)  and  the  model  by 


Sodini  et  al.  (1984),  are  quite  similar.  This  may  mean  that  there  could  be  some 
uncertainly  in  the  device  parameters  given  by  Chung  et  al.  (1991). 


vr.r-  (V--) 


Fig.  6.  Plot  for  the  determination  of  the  parameters  p  and  Vl. 

An  alternative  way  to  crudely  estimate  p  from  the  experimental  data  (often  used 
in  literature)  is  to  consider  the  saturation  current  slightly  above  threshold  where 
(VgtI^l)^  «  1-  Then,  the  equation  for  the  saturation  current  can  be 
approximated  by 


and  (p/2)^/^  is  found  as  the  slope  of  the  plot  of  Usat  versus  Vgt  (=  ^gt  in 
this  regime).  (It  may  also  suffice  to  measure  Isat  for  one  value  of  Vgt  sufficiently 
close  to  threshold).  However,  this  technique  is  not  very  accurate  in  short  channel 
devices  where  the  region  Vgt  «  Vl  is  very  narrow  and  may  be  affected  by 
subthreshold  current  caused  by  short  channel  effects. 

Using  the  various  techniques  described  above,  we  arrive  at  the  following 
parameter  values  used  to  fit  the  experimental  l-V  characteristics  in  Figs.  1  and  2: 
3  »  0.0067  and  Vl  *  0.36  V,  which  correspond  to  *  4x10^*  m/s  and 
|i/i  *  0.028  m^/Vs.  However,  we  note  that  a  higher  value  of  3  (and  pn)»  i  e., 
more  in  accordance  with  the  values  obtained  from  the  slope  in  Fig.  4,  would  have 


given  a  better  fit  in  the  linear  part  of  the  l-V  characteristics,  while  the  present 
values  favour  a  good  fit  in  the  saturation  region  (which  is  of  primary  importance 
in  digital  applications). 

Finally,  we  have  to  choose  the  voltages  V(^(  and  Vc  which  determine  the 

dependence  of  the  threshold  voltage  parameter  o  on  the  gate-source  voltage. 
Here,  we  have  chosen  Vc  =  0.2  V  and  Vac=  Vjo  +  0.5iVgi)max  =  1-7  V. 

However,  these  parameter  values  can  be  slightly  adjusted  for  a  better  fit. 

4.  Application  of  Model  in  AIM-Spice  Simulation 

The  present  model  has  been  implemented  in  our  circuit  simulator. 
Automatic  Integrated  Circuit  Modeling  Spice  (AIM-Spice).  AIM-Spice  runs 
under  the  WINDOWStm-3  environment  on  IBM  PCs  and  compatible  computers 
and  takes  full  advantage  of  the  resident  graphics  user  interface  and  memory 
above  the  640k  DOS  limit. 

Fig.  7  is  a  screen-dump  from  the  AIM-Spice  post-processor  showing  a 
composite  view  of  relevant  voltages  and  currents  from  a  turn-off  transient 
analysis  of  a  simple  NMOS  inverter  circuit  with  an  NMOS  load  at  the  output. 
The  subthreshold  regime  is  emphasized  by  showing  the  turn-off  recovery  of  the 
drain  current  of  the  driver  transistor  in  a  semi-logarithmic  scale.  The  inverter 
and  load  contain  altogether  three  identical  quarter  micron  transistors  with  the 
specifications  corresponding  to  the  transistors  with  characteristics  shown  in  Fig. 
1.  The  capacitances  are  calculated  according  to  Meyer's  expressions  (see  Section 
2). 


Fig.  7.  Transient  analysis  of  an  NMOS  inverter  with  an  NMOS  load  at  the  output.  All  transistors 
have  the  specifications  given  in  Fig.  I .  The  trace  in  the  background  shows  the  input  to  the  gate  of 
the  driver  transistor  (the  gate  of  the  inverter  load  transistor  is  fixed  av  15  V  and  the  voltage  supply 
is  Vdq  =  3.0  V).  The  upper  and  lower  traces  on  the  right  hand  side  show  the  drain-source  voltage 
and  the  drain  current  for  the  driver  transistor,  respectively,  (from  Ytterdal  et  al.  ( 1991): 


III.  Modeling  of  Radiation  Effects. 


Fig.  8  shows  typical  experimental  data  for  an  imhardened  CMOS  process 
(from  Baze  et  aL  (1990)).  As  can  be  seen  from  the  figure,  the  60  krad  dose  leads 
to  a  considerable  shift  of  the  threshold  voltage  and  to  a  drastic  (over  three  orders 
of  magnitude)  increase  in  the  leakage  current.  The  device  leakage  current 
recovers  after  a  long  aimeal.  However,  the  threshold  voltage  shift  remains. 


Fig.  8.  Current-voltage  characteristics  for  unhardened  CMOS  (cfter  Baze  et  al.  (1990)). 

n-channel  device,  Vgs-5V,  Dose  rate  167  rad/s.  0  -  pretest,  2-60  krad  (Si02), 
1  60  krad  and  107  sec  anneal  at  25  ®C. 

These  characteristics  can  be  readily  modelled  using  our  model  when  we  include 
tm  additional  leakage  current 


^leak  —  (Sieak^ds  (32) 

The  calculated  curves  are  shown  in  Fig.  9. 
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Fig.  9.  Calculated  current-voltage  characteristics  for  unhardened  CMOS.  0..  ■  pretest,  I .- 60 
kr^  after  anneal,  2. -60  krad  dose 


IV.  Software  listing. 


REM 

REM  Michael  Shur  (804)-924-6109 
REM  fax  (804)-924-88l8 
REM  e-mail  ms8n@virginia.edu 
REM 

REM  Program  MOSFET.s 
REM  version  I  5/2/91 

REM  pMOSFET.s  calculates  unified  short-channel  MOSFET  I-Vs 
REM 

REM  The  program  is  the  subroutine  of  program  PLOTII 

REM  which  creates  data  files  for  Cricket  graph 

REM 

REM 


SetLimits: 

REM  x-axis  variable  limits 
Xstart  =-l:Xend=2 
RETURN 
ReadData; 

Fpar%=l 

A$="pMOSFET.s" 

FL$="pMOSFET.s'’ 

REM  +++++  DATA  SECTION  ++++++ 

REM  constants 

REM  -H-+++  DATA  SECTION  i 

REM  constants 

q=1.602E-19  Coulomb 

kb=1.38E-23’J/m 

REM  temperature 

T=300'  degree  K 

REM  material  parameters 

eps=1.04E-10'F/m 

epi  =  3.44E-irF/m 

W  =  .00003 

L  =  .0000007 

di  =  .0000001  . 

Nd  =  4E+23 
mob  =  .028 
'm2/Vs  mobility 

vs  =  40000!’m/s  saturation  velocity 

m=4'  saturation  knee  parameter 

Vto  =  .44  'threshold  voltage 

Rsm  =  .75  'source  series  resist,  ohm-mm 

Rdm  =  .75  'drain  series  resist,  ohm-mm 

Rs=Rsm/(W*1000) 

Rd=:Rdm/(W*1000) 

Rt  =  Rs  +  Rd 

Iambda=.044'  output  conductance  parameter 

sigmao=.048  'DIBL  parameter 

Vsigma  =  .2  'DIBL  voltage  parameter 

Vst  =5  1.75  'DIBL  transition  parameter 

eta=1.32  'subthreshold  slope 

sig=lE-ll 

Vds=l 

GOSUB  SetLimits 
PRINT  A$:  PRINT 


GOSUB  PARDEF 

RETURN 

PARDEF: 

Jpar%(l)  =  3:par(l,l)=0:par(  l,2)=l;pur(1.3)=2 

Rad  =  PARAMETER 

RETURN 

UPDATEPAR: 

GOSUB  PARDEF 

RETURN 

FirstFunction: 

Vgs  =  XXX 

IF  Rad=0  THEN  Vto=.5;sig=8E-l2:eia=l  .72 
IF  Rad=l  THEN  Vto=.25:sig=I.2E-l  l:eta=L72 
IF  Rad=2  THEN  Vto=.l  :sig=.0000()()l  :eta=6 
Vth  =  kb*T/q'  thermal  voltage 
VL  =  vs*L/mob 
no  =  epi*eia*Vth/(2*di*q) 

Vgto  =  Vgs  -  Vto 

sigma  =  sigmaoA  l+EXP((Vgto-Vst)/Vsigma)) 

Vgt  =  Vgto  +  sigma*Vds 
Vgho  =  .5*Vgt/Vth 
Vgh=:Vgt/(eta*Vth) 

Vgte=Vih*(  I +Vgho+(9+(  Vgho- 1  )''2)''.5 ) 

GOSUB  nscalc: 

'PRINT  Vgs,  Vgte,  Vgh.  ns 
gchi  =q*ns*W*mob/L 
gch  =  gchi/(l+gchi*Rt) 
gR=gchi*Rs 

Isat  =  gchi*Vgte/(l-hgR+SQR(l-*-2*gR+(Vgte/VL)''2)) 
zz=gch*Vds/Isat 

Ids=gch*(I+Iambda*Vds)*Vds/(I+zz''m)''(l/m)+sig*Vds 
yyy  =  LOG(Ids)/LOG(IO) 

PRINT  sig,  Vto,eta 
RETURN 

REM  subroutine  nscalc 
nscalc: 

REM  ns  calculation 

IF  Vgh>  86  THEN  ns  =  epi*Vgt/di 

IF  Vgh<  -12  THEN  ns  =  no*EXP(Vgh) 

IF  Vgh>  -12  AND  Vgh<86  THEN  ns  =  2*no*LOG(l  +  .5*EXP(Vgh)) 
RETURN 
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VII.  Appendix.  Device  structure  and  equivalent  circuit 


Fig.  Al.  Schematic  device  structure  of  NMOS. 

Typical  device  dimensions:  W  =  30  pm,  L  =  0,7  pm,  di  =  100  A 
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Fig.  A2.  Sinq)lified  equivalent  circuit  of  NMOS. 


